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Abstract 

We use a hydrodynamic model to describe the relaxation of optically injected currents in quantum 
wells on a picosecond time scale, numerically solving the continuity and velocity evolution equations 
with the Hermite-Gaussian functions employed as a basis. The interplay of the long-range Coulomb 
forces and nonlinearity in the equations of motion leads to rather complex patterns of the calculated 
charge and current densities. We find that the time dependence of even the first moment of the 
electron density is sensitive to this complex evolution. 
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The dynamics of hot electron currents, which determines the distance that injected car- 
riers can propagate and the charge density patterns that can form, is important for the 
design of various semiconductor devices, including transistors, charge-coupled light detec- 
tors, cascade lasers and light emitting diodes 1 . The full dynamics of systems strongly out 
of equilibrium requires a complicated analysis based either on extended numerical modeling 
employing Monte Carlo simulations 2 , or on quantum versions of the kinetic equations^ 1 ^^. 
For two-dimensional (2D), systems quantum kinetic approaches are numerically demanding 
and often unfeasible. 

The direct injection of strong electrical and spin currents^ is possible using the interference 
of one-photon (frequency 2uo) and two-photon (frequency u) absorption. The high speed of 
the injected electrons, v ~ 1000 km/s, is determined by the excess photon energy 2hu — E g , 
where E g is the band gap for the quantum well (QW). In an experiment injecting electrical 
current, coherent control is achieved by adjusting the relative phase parameter A<j) = 20 w — 
4>2u) of the two beams, where W (02^) is the phase of the beam at u> (2a;) . The resulting 
current density is approximately given by Nv sin(A0), where N is the concentration of 
optically injected carriers. 

In experiments done on multiple quantum well (MQW) samples, with the beams incident 
along the growth direction, these injected lateral currents (see Fig.l) have been detected 
using various techniques^. The QWs are well-separated and unbiased, so electrons do not 
have sufficient energy either to tunnel through the barriers or overcome them by thermal 
activation. For this reason no dynamics along the "vertical" growth direction is possible, and 
in these systems we do not explore the rich nonlinear behavior seen in vertical transport in 
superlattices.— . Instead, our subject is only the lateral dynamics of the carriers subsequent 
to their injection. 

Nonetheless, the analysis of this lateral dynamics is more complicated than in many 
transport problems because of the largely inhomogeneous long-range Coulomb field, leading 
to complicated space-charge effects. As the electrons and holes move away from each other, 
the characteristic range of inhomogeneity is on the order of the size of the electron and 
hole puddles themselves. A qualitative approach to this problem is to adopt a "rigid spot" 
approximation, which reduces the problem to the motion of two coupled and damped linear 
oscillators, representing the centers of the electron and hole puddles, each moving with 
a uniform velocity and thus exactly preserving its shape.— Here the entire dynamics is 
described by four parameters: displacements of electron and hole spots and their velocities. 

An important question is to what extent this "rigid spot" model can describe the dynam- 
ics, as Coulomb and other concentration-dependent forces come into play when the puddles 
separate. Typically, the change in carrier densities at any point is small, since the relaxation 
and Coulomb interaction prevent any large separation of the centers of charge. However, 
can the spatial dependence of this small change depart significantly from what is predicted 
in the "rigid spot" model, leading to a picture involving more generally "distorted puddles" 
as more appropriate ? 
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Answering this question is important for the ultimate design of numerical schemes to 
efficiently explore the system dynamics. What is needed initially is a practical model and 
calculation procedure that captures the essential physics, is transparent enough to allow 
extraction of the important features of the dynamics, and reasonable enough to allow com- 
parison with experimental results, at least at a semiquantitative level. Here we propose a 
phenomenological hydrodynamic model for the dynamics of optically injected charge currents 
subject to space-charge effects. The advantage of the hydrodynamic approach is that it is 
insensitive to the details of the carriers' distribution functions, with the dynamics described 
by a coupled system of partial differential equations for the concentrations and velocities. 
The collisions between carriers, and interaction between the carriers and the environment, 
are modeled by a set of characteristic times determined by the experimental conditions. 

As an example calculation we consider a MQW structure (Fig. 1), assuming the excitation 
and subsequent dynamics in each quantum well is the same; within our model we perform a 
full 2D calculation of the carriers' motion. The u- and 2a;-beams produce initial distribution 
of electron and hole densities and velocities which then evolve in time and space. The initial 
distribution of carriers in each single QW is Njfl (r, t = 0) = N s exp (— r 2 /2A 2 ), where A 
is governed by the beam sizes (see caption of Fig.l). The total concentration N e<h (r, t) = 
qNffl (r,t), where r = (x,y) is the 2D coordinate, and q is the total number of single QWs, 
neglecting the distance from the first to the last single QW in the structure compared to A. 
The spin indices are omitted since here the electrons are assumed to be unpolarized. The 
(r, t) arguments will be omitted below for brevity. 

For simplicity, in this letter we present results for the dynamics of electrons assuming 
that the holes are infinitely heavy; we will include hole dynamics in a later publication. The 
equations describing the motion of carriers consist of the continuity, and momentum and 
energy evolution sets. In the effective mass approximation, the analysis of the dynamics 
shows that a first description can be provided even without taking the energy relaxation 
into account. For this reason we restrict ourselves to the first two sets of variables. The 
continuity and the evolution of the local mean velocity u (r, t) of electrons, are described by: 



where E = E e + E is the macroscopic Coulomb field produced by electrons and holes, and 
e is the elementary charge. Here m e is the electron effective mass, r e h describes momentum- 
conserving drag due to the Coulomb forces at electron-hole collisions^ and can be weakly 
concentration-dependent itself^, r e is the relaxation time for electrons due to external fac- 
tors, such as phonons and impurities leading to the relaxation of the total momentum, 
and Nq = qN s . In addition to the explicit time scales r e and r e h a third important time 
scale in the problem is where Q p i = (7r/2) 3//4 J N e 2 / em e ^/2A is the characteristic two- 
dimensional plasma frequency^ (e is the dielectric constant); it characterizes the strength 
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Figure 1: (Color online) A sketch of the optical current injection scenario. The MQW structure 
with q = 4 (gray bands for GaAs) is shown. Vectors eu2u>) correspond to the polarization of the 
photons. The injected charge density follows the profiles of the 2a;-beam and the square of the 
w-beam intensity, which are the same due to the 2tu generation procedure (frequency doubling in 
nonlinear crystal). U(0) is the initial macroscopic speed of the electron puddle. 

of the Coulomb interaction. In our sample calculations we take r e = 80 fs and r e h = 150 fs; 
a typical value N s = 10 11 cm -2 with q = 1 and our choice of A = 1 fxm leads to Qi~f « 1.4 
ps. These times agree with the set of parameters of Due et al. [4| and are shorter than those 
for vertical transport 1 ^ due to the higher dimensionality of carrier motion here. 

To reduce the system of partial differential equations, we introduce a full basis set in the 
form of the harmonic oscillator wave functions 

(r) = VW(*M»(V), M*) = ^77T= e ' x2/2A ^ ( 3 ) 

where H n (x/A) is the Hermite polynomial of the n— th order. The key point of our approach 
is the expansion of the possible solutions in a finite basis in the form: 

N e = E ^i^(*)*n ltn9 W ; U i= J2 to + U { (t) . (4) 

ni,n,2 ni,ri2 

Here i = x, y is the Cartesian index. To improve the convergence, we include a known 
function of time U{ (t) in the right-hand side of Eq.(4). This function is obtained by solving 
the linear equations of motion in the rigid spot approximation^. In the geometry considered 
here, U y (t) = 0, and, therefore, we drop the Cartesian index of U x (£). The electric fields at 
the point r are 

(r,t) = T- E / ( ? ) r^s^ ( 5 ) 

£ — J r — r 

ni,n 2 I I 

where the upper (lower) sign corresponds to electrons (holes). If a QW is located close 
to the semiconductor-air interface, the role of the image charges which can be taken into 
account with the Green function technique 14 , increases the field by a factor of two. Since 
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A is typically much larger than the distance between the QWs, within our assumption that 
the excitation of all the wells is the same, the electric field at a point r in each well and the 
subsequent dynamics is the same. 

The state and the motion of the electron spot is then fully described by a 3 x 
("-max + l) 2 —component vector S a , where a is a two-component index corresponding to 
indices n x and n 2 . By projecting the equations (1),(2) and (5) on the set of ^ ni ,n 2 ( r ); we 
obtain the system of ordinary nonlinear differential equations in the form 

-Jp = ^2C(a;r],fx)S ri S IM - ^r(a;C)S c . (6) 

w C 

Here the matrix C(a;r],n) is determined by the spatial dependence of concentration and 
velocity, and the electron-hole drag, while the matrix T (a; Q depends on the Coulomb 
forces and momentum relaxation times. The full matrices will be presented elsewhere. 

The initial conditions are given by: (i) Aq (0) = y/nN , N* n (0) = (if n\ > or 
n 2 > 0), and (ii) Unl% 2 (0) = 0. Conditions (i) correspond to the initial injection of density 
in the Gaussian mode \I/ ,o ( r ) only. Due to the symmetry of the problem 2m+1 (t), 
u ni2m+i(t) an d u ni,2m(t) wm remain zero. The initial speed of the electron spot U(0) is 
determined 7 by the photon excess energy and A(f>. 

At very short times after the injection the carriers propagate ballistically, and then the 
motion becomes diffusive and influenced by the Coulomb forces. A simple characterization 
of the motion of the charge is given by the first moment of the concentration, 

(x(t)} = ±-JxN e (r, t) d 2 r, N t = 7rA A 2 , (7) 

where N t is the total number of injected electrons. The relaxation and Coulomb interaction 
guarantee that £ = (x(t)) /A is always small in our systems. Were the rigid spot approx- 
imation valid, the amplitude of the nonzero higher terms N£ Q would drop off as £ n . The 
main qualitative result of our investigations is that this is often not so for reasonable exci- 
tation parameters. Instead, the concentration-dependent Coulomb forces and drag lead to 
a deformation of the spot and, therefore, to an increase in the number of ^ ni ,n 2 terms with 
non-negligible amplitude in the sums, and to a complicated spot structure. In the course of 
time, a few of the n2 (t) components can become of the same order of magnitude, although 
all small compared to A 0j o(0) due to a relatively small spot displacement ^< 1. The number 
of higher terms that are important depends mainly on the structure of the C (a;rj, /^-matrix, 
with their relative contribution dependent on the parameters T e Q p \ and T e /r e h. In contrast 
to the scattering by phonons and disorder, which stabilizes the spot motion and does not 
lead to its distortion, the role of electron-hole collisions is two-fold. On one hand, it stops 
the motion of the spot as a whole, on the other hand, it generates higher components in the 
velocity pattern due to nonuniform velocity evolution. The nonuniform velocity distribution 
produced by the Coulomb forces and the drag then leads to a nonuniform density distri- 
bution due to continuity. So the joint effects of the Coulomb forces and the drag lead to 
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Figure 2: (Color online) The change in the electron density distribution compared to the injected 
density at t = 80 fs (left column) and t = 1.5 ps (right column). The MQW parameters are q = 1 
(upper row) and q = 16 (second row). Relaxation times are r e = 80 fs, T e h = 150 fs. The spot size 
A = 1 fim, N s =10 n cm" 2 , 17(0) = 400 km/s. 

a "buckling" of the electron spot along the direction of its velocity. The buckling becomes 
more pronounced as the space-charge effects, characterized by the parameter T e Q p i, increase. 

This is illustrated in Fig. 2, where we show the pattern of the change in the electron density 
at two times, for samples with q = 1 and q = 16 buried deep within bulk Al x Gai„ x As where 
Eq.(5) is applicable. On the time scales of the ballistic regime, t < r e , in both samples the 
spot moves as a rigid distribution with 

N e = N exp [- (x - (x(t)}) 2 /2A 2 ] exp [-y 2 /2A 2 ] , (8) 

and the corresponding density component develops as Nf fi (t) = A^ (0)^/v^2- This is clearly 
seen in the left column of the Fig. 2, where the density profile is almost the same for systems 
with weak (q = 1) and much stronger (q = 16) space-charge effects at t = 80 fs, and is well 
described by Eq.(8). At 1.5 ps, the profile for the q = 1 system buckles only slightly; a 
"rigid spot" picture is applicable. But for q = 16 the behavior at these times is much 
more complicated, with several ^ ni ,n 2 states contributing, and the profile changing from 
oval to bow-shaped. Here the "rigid spot" approximation is no longer reasonable, and a 
more general "distorted puddle" picture emerges. More insight can be obtained from the 
current density N e (r)u x (r) distribution in Fig.3. At t = 80 fs the patterns for q = 1 and 
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Figure 3: (Color online) The electron current density distribution (in arb. units) at t = 80 fs (left 
column) and t = 160 fs (right column). The MQW parameters are q = 1 (upper row) and q = 16 
(second row). Other parameters are the same as in Fig. 2 

q = 16 look similar. At t = 160 fs, the distributions already differ very considerably, and 
the pattern for q = 16 is much more nonuniform than that for q = 1. 

It is interesting to see the effect of the real spot shape on gross features of the density 
distribution, such as its first moment (j7|). We evaluate this from our calculations as 

^lim ^max JSJ^ (j^^) f f 

( x (t)) = Yl Yl ~^/\r — / x ^ ni ^ dx / ^ m ^ dy ' ( 9 ) 

ni=ln 2 =0 * J J 

and present in Fig. 4 the results for q = 4 and q = 16 samples. The upper limit ri\ im here is 
less than the number of Hermite- Gaussian states in the basis n max , since the contributions 
of the upper states cannot be calculated with a high precision; this problem is typical for 
calculations with finite basis.— For this set of parameters we find that n max = 14 and 
n Um — 7 give converged results. In addition, Fig. 4 shows the displacement in the rigid 
spot approximation^, where the dynamics yields (x Tsa (t)) = £/(0)e~ 7 *sinh (7 P it) /j v i with 
27 = l/r e + l/2r e /i, 7^ = 7 2 — fipj. With the stronger space-charge effects that arise for 
a larger number of wells, the deviation from the rigid spot approximation becomes more 
severe. However, both (x TSa (t)) and (x(t)) reach their maxima at almost the same time 
to = 2r e |ln(O p iT e )|, and the maximum values are very close to each other. The difference 
begins developing at t > to, when the buckling effects become more pronounced. The curves 
for small value of r e h = 50 fs show that a decrease in r e h does not qualitatively modify the 
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time [ps] 

Figure 4: (Color online) Mean spot displacement (x(t)) (solid lines) and (x TSSb (t)} (dashed lines). 
Parameters q and r e h are marked near the lines. Other parameters are the same as in Fig. 2. 

dynamics in our range of parameters. However, in a regime with r e h <t e a decrease in r e h 
can lead to an increase in the puddle distortion. We turn to the influence of this distortion 
on experimental results in a separate publication. 

To conclude, we have shown that within a hydrodynamic model a Hermite-Gaussian 
expansion can be used to study the dynamics of currents injected optically by coherent 
control in MQW samples. As samples with increasing numbers of QWs are considered, we 
find the resulting dynamics shows a complicated buckling behavior in the density distribution 
that is driven by space-charge effects, with the charge density departing more and more from 
what a rigid spot approximation would predict. The details of our calculated results reflect 
both the simple dynamics (1),(2) we adopt, and the assumption of equal excitation of all the 
quantum wells. But it is unlikely that a relaxation of these approximations, and the richer 
dynamics that would result, would lead to a more rigid behavior of the charge puddle. 

For T e h = 150 fs, we find that the complex behavior, associated with a number of com- 
parable Hermite-Gaussian amplitudes, emerges at fl p \T e « 0.1, well within the overdamped 
regime, and is determined by this parameter rather than by the electron-hole separation. 
How sensitive various experimental probes will be to this deviation from rigid spot dynamics 
has yet to be examined; at least it has significant consequences on the evolution of the first 
moment of density on the timescale of a picosecond. 

Our numerical consideration relies on the finite basis analysis. In the overdamped regime 
we have considered here a fast decay of the off-diagonal Coulomb matrix elements with 
the distance between the states ensures the applicability of relatively small basis sets. As 
we consider systems with stronger Coulomb interaction, the basis should be extended. We 
can naturally expect the distortion of the carrier puddles to get more pronounced, and an 
interesting problem is whether it can lead to a chaotic behavior— when many Hermite- 
Gaussian terms form the density and velocity patterns. 
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